Surface tension of electrolytes: Hydrophilic and hydrophobic ions near an interface 
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We calculate the ion distributions around an interface in fluid mixtures of highly polar and less 
polar fluids (water and oil) for two and three ion species. We take into account the solvation and 
image interactions between ions and solvent. We show that hydrophilic and hydrophobic ions tend 
to undergo a microphase separation at an interface, giving rise to an enlarged electric double layer. 
We also derive a general expression for the surface tension of electrolyte systems, which contains 
a negative electrostatic contribution proportional to the square root of the bulk salt density. The 
amplitude of this square-root term is small for hydrophilic ion pairs, but is much increased for 
hydrophilic and hydrophobic ion pairs. For three ion species including hydrophilic and hydrophobic 
ions, we calculate the ion distributions to explain those obtained by x-ray reflectivity measurements. 

PACS numbers: 



I. INTRODUCTION 

It has long been known that the surface tension 7 of 
a water-air interface increases with addition of inorganic 
salts in water. Wagner— ascribed its origin to the image 
charge in air due to the difference in the dielectric con- 
stants of the two phases, which repels each ion in water 
away from the interface. Onsager and Samaras^ obtained 
the limiting law for the surface tension change in the form 
A7 = ■|Tn„^B[ln(l/nw^B)+const.], where l^^) is 

the ion density in the bulk water and is the Bjerrum 
length. Levin and Flores-MenEfc^ took account of ion de- 
pletion due to a relatively large size of the hydration shell 
radius^. In experiments, at not extreme dilution, the 
linear behavior A7 = TnwAg has been measured^i^iii^, 
where is the effective thickness of the ion-free layer. 
For example, A^. ^ SA for NaCl, where the densities of 
Na"*" and Cl^ are n^/S. However, for a number of salts 
around 1 mM in aqueous solutions, Jones and Ray^ de- 
tected a very small negative minimum in A7, which still 
remains an unsolved controvers y^' ^ ' . 

In polar fluid mixtures, the ion distributions are much 
more complicated when the ions are soluble in the two 
phasesiiii^iii^. In this problem we need to account for 
the solvation between ions and solvent molecules (hy- 
dration in aqueous solutions), whose free energy contri- 
bution usually much exceeds the thermal energy T (per 
ion)^. In particular, one ion species can prefer one fluid 
component, while the other species can prefer the other 
component. Such asymmetric ion pairs tend to segre- 
gate around an interface, while the segregation is pro- 
hibited on larger scales due to the charge neutrality in 
the bulk. In this paper, we will examine this microphase 
separation near an interface to calculate a decrease of 
the surface tension A7 < 0, whose amplitude can be 
much larger than that of the well-known increase for 
hydrophilic ion pairs. It is worth noting that, in their 
small- angle neutron scattering experiment, Sadakane et 
ali^ found periodic charge-density-wave structures in 
a near-critical binary mixture of D20-trimethylpyridine 
containing sodium tetrarphenylborate (NaBPh4). Their 



salt is composed of strongly hydrophilic cation Na+ and 
strongly hydrophobic anion BPh4^, which should consid- 
erably decrease the surface tension and produce a meso- 
scopic structure near the critical point. 

In their x-ray reflectivity experiment, Luo et al^ mea- 
sured the ion distributions in the vicinity of an interface 
in water-nitrobenzene. They added two salts, tetrabuty- 
lammonium tetraphenylborate (TBA-TPB) and tetra- 
butylammonium bromide (TBA-Br). Then they realized 
a two-phase state with hydrophilic anion Br~ , hydropho- 
bic cation TBA+, and hydrophobic anion TPB^. They 
detected strong accumulation or depletion of the ions on 
the two sides of the interface, which suggest a crucial role 
of the ion-solvent interactions dependent on the ion and 
solvent species. In water-nitrobenzene containing both 
hydrophilic and hydrophobic ions, a large drop of the 
surface tension has been observe d^^i^^ . 

We mention the presence of a large body of spec- 
troscopic experiments and computer simulations with 
molecular resolution on ion effects at a water- air 
interfacei^iii, where the air region may be treated as 
a vacuum. Such studies provide detailed information 
of ionic interfacial behavior on the angstrom scale for 
various ion species. However, microscopic studies re- 
main inadequate for ion effects at a water-oil interface. 
On the contrary, our approach in this paper will be 
based on a Ginzburg-Landau theory of solvation and ion 
distributionsi^ii^. We will consider water-oil systems like 
water-nitrobenzene, where the dielectric constants of the 
two components are not much separated (the dielectric 
constnt of nitrobenzene is about 35). 

The organization of this paper is as folows. In Sec- 
tion 2, we will present a Ginzburg-Landau approach 
to the molecular interactions between ions and solvent 
molecules. Taking account of the electrostatic, solvation, 
and image interactions, we will introduce the grand po- 
tential and present a theoretical expression for the surface 
tension. It contains a negative electrostatic correction, 
which is not included in the Gibbs theor y^^i^^ but is cru- 
cial for hydrophilic and hydrophobic ion pairs. In Section 
3, we will numerically examine the ion distributions and 



2 



the surface tension for hydrophihc and hydrophobic ion 
pairs. We wiU also discuss the behavior of A7 in the 
usual case of hydrophihc ion pairs in extreme dilution. 
In Section 4, we extend our theory in the presence of 
three ion species to explain the experiment by Luo et 
al^. In Appendix A, we will examine the relationship 
of our surface tension formula and the Gibbs equation. 
In Appendix B, we will relate our grand potential to the 
bulk pressure and derive the Laplace law using our sur- 
face tension formula, since we will treat incompressible 
binary mixtures in the text. 



II. THEORETICAL BACKGROUND 

A. Ginzburg-Landau free energy 

We consider a polar binary mixture containing a small 
amount of salt. The volume fraction of the more po- 
lar component is written as (p. The other less polar 
component has the volume fraction 1 — (p. We neglect 
the volume fractions of the ions. The ion densities are 
written as ni,n2, • • ■. Their charges are Z^e, so .Zi = 1 
and Z2 = — 1 for two monovalent ion species. In our 
scheme, (p, ni,n2, • • • are smooth space-dependent vari- 
ables coarse-grained on the microscopic level. We ex- 
amine the ion distribution around an interface^, where 
all the quantities change along the z axis. The sol- 
vent molecular sizes are given by a common length a 
for the two components. Then (p also represents the mo- 
lar concentration. The interface thickness ^ is assumed 
to be longer than a. We neglect the formation of dipole 
pairs and ion clusters, which is relevant at not small ion 
densitiea^SiSiiS^. 

The free energy F of our system is the space integral 
of the free energy density / of the form, 



/ = f,{cp,T) + ^\\/(t>\-' + 
+T ^ Ui [In(nia^) - gn 



g(0) 



(2.1) 



We set the Boltzmann constant equal to unity. The first 
term /o is the chemical part in the Bragg- Williams form. 



/o = ^[01n(/.+ (l-<^)ln(l 



+ XHI -<!>)], (2.2) 



where x depends on the temperature T and its critical 
value is 2 in the absence of ions^^'^'*. The second term is 
the gradient part, while the third term is the electrostatic 
free energy with (f> being the electric potential. The elec- 
trostatic potential $ satisfies V • e{(j))W^ — — 47r/3, where 
p = Zierii is the charge density. Around an interface 
the electric field E = —d^/dz is expressed as 



E{z) 



47r 



e(0(z)) 



dz'p{z'), 



(2.3) 



where the lower bound of the integration is pushed to 
—00. The dielectric constant e is assumed to be of the 
linear form^^ 



e{(p) =£c + £i{(p - (pc), 



(2.4) 



where Sc and Si are constants. The cpc is the crit- 
ical volume fraction (= 1/2 for the free energy den- 
sity in Eq.(2.2)). Thus e{(j>) depends on z near an 
interface. The coupling terms —TgiUiCp arise from 
the composition-dependence of the solvation chemical 
potentials'^. The differences Tg.iA<p may be equated to 
the Gibbs transfer energies (per particle here) known 
in electrochcmistryiii^i^li^iS^, where Acj) is the con- 
centration difference between the two phases. From 
data of room-temperature water-nitrobenzene in strong 
segregation'^-, gi are expected to be typically of order 15 
for monovalent hydrophihc ions and are even larger for 
multivalent ions such as Ca^"*" or Al^^^ . For the hydropho- 
bic anion BPh^ (tetraphenylborate) , for example, it is 
about —15, on the other hand. The resultant solvation 
coupling between the ions and the composition is thus 
very strong, dramatically affecting the phase transition 
behavior near the critical pointi^ii^. In binary mixtures 
in which the more polar component is water, the cou- 
pling constants gi are positive for hydrophihc ions and 
negative for hydrophobic ions. 

The last term in Eq.(2.1) represents the image inter- 
action, arising from inhomogeneous dielectric constant 
e. It originates from the discrete nature of ions, while 
the electric field $ in our theory is produced by the 
smoothly coarse-grained charge density p. The origi- 
nal papersi^ treated water-air interfaces, but we sup- 
pose weak or moderate inhomogeneity of the dielectric 
constant across a diffuse interface. Then it follows the 
Cauchy integral fornii^. 



1k\z — z'\ 
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d(p{z') 



dz' 



(2.5) 



to first order in ei. The coefficient A represents the 
charge strength as 



(2.6) 



where £bc — e'^/^cT is the Bjerrum length at e = Ec- 
The damping factor g-'^'^l^-^ I in Eq.(2.5) arises from the 
screening of the potential by the other ionaii^. In our 
numerical analysis, we treat n as the space-dependent 
local value [47re^m(r)/e(r)r]'/^ with m ~ Zfui. 

In equilibrium we assume homogeneity of the chemical 
potentials pi = SF/Srii {i = 1,2,---) and h = 6F/6(p. 
We introduce normalized ion densities^, 
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a n, 



(z = l,2,...). 



(2.7) 



In our theory the ion volume fractions are assumed to 
be so small such that the real ion sizes do not come into 
play. Since we have 



fit = T(lncj + 1- g^' 



y2 



(2.8) 
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the ion density profiles are expressed as 



exp 



Zje 



T 



A*im(2) 



(2.9) 



where = exp(/ii/r — 1). In taking the derivatives 
of the image free energy — J dr^i^ J2i ^f'f^i with 
respect to rii, we neglect the n-dependence of /iim in 
Eq.(2.5). The h = SF/S4> has the meaning of the chem- 
ical potential difference of the two fluid components (di- 
vided by a^). It is of the form, 



fill 



(2.10) 



where /q = dfo/dcj) and 0" = cPip/dz'^. From the image 
interaction we have the contribution /lim — SFi^/Scf). In 
the ID case, /ii,„ is given by the right hand side of Eq.(2.5) 
with dtp{z')/dz' being replaced by Zfdni{z')/dz'. 

Hereafter the quantities with the subscript a (the sub- 
script (3) denote the bulk values in the more (less) polar 
phase attained as z — )■ —oo (as z oo). We assume that 
the screening lengths, and much shorter 

than the system length L. In real systems this might not 
be the case for extremely small np. From Eq.(2.9) the 
distribution ratio of each ion can be expressed as 



Ctp/cia = exp[ZjeA$/r - g^A^], 



(2.11) 



where A0 = 4>a — 'Pp- The image interaction vanishes 
far from the interface and does not appear in Eq.(2.11). 
There arises a Galvani potential difference A$ = $0 — 
across an interface. It is determined by the charge 
neutrality far from the interface, 



(2.12) 



In addition, as z ±cxd, the homogeneity of h in 
Eq.(2.10) yields the bulk relations. 
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(2.13) 



-^shein (o'^ order of the size of a water molecule). In 
our scheme, for finite interface thickness ^, hydrophilic 
ions are repelled from an interface into the a region due 
to the solvation interaction (due to the factor e^'"^ in 
Eq.(2.9) for gi > 0). Thus, even in the absence of the 
image interaction, a depletion layer of hydrophilic ions 
can be formed and the linear behavior A7 cx still 
follows. To make qualitative arguments, therefore, the 
image interaction may be neglected for not very large A 
in Eq.(2.6). 




FIG. 1: Image factor -Fim(z) for monovalent ions around an 
interface for (a) x = 2.3 and A — 4, (b) x = 3 and A — 4, 

and (c) X = 3 and A = 10, where cia = cifs = 2 x 10"''. 

In Fig. 1, we show examples of the image factor for 
monovalent ions near an interface, 



Fhn{z) = exp[-^i,„(z)/r]. 



(2.15) 



which is calculated from Eq.(2.5) for Ei/Ec = 0.8 and 
X = 3. This factor appears in Eq.(2.9) for Zi — ±1. It is 
smaller than unity in the a region and is larger than unity 
in the f3 region, since the image potential is repulsive 
in the a region and is attractive in the /3 region (see 
Appendix B of our previous paperA'^). At low densities, 
its minimum and maximum are more enhanced for larger 
A ~ iBc/a and for larger x (larger At/)), in accord with 
Eq.(2.14). See also analysis of the image interaction in 
Figs. 8-10 of our previous worki^. 



B. Image interaction 

We give more discussions on the image interaction. 
From Eq.(2.9) it is important under the conditioni^, 

^ < z2£Bc£iA0/4e, < (2.14) 

where k = k„ or Kp in the two phases. When ions are 
present only in the more polar phase under Eq.(2.14), the 
change of the surface tension obeys the Onsager-Samaras 
law^. In the thin interface limit ^ ^ 0, Levin and Flores- 
Mena'^ argued that ions in water cannot approach the 
interface within the distance of the hydration shell radius 



C. Surface tension 

In order to calculate the surface tension, we introduce 
the grand potential density, 

uj = f -'^^liUi- h(j), (2.16) 

i 

where / is given by Eq.(2.1). The space integral of cj 
is minimized under given boundary conditions in equi- 
librium. In Appendix B, we will examine how uj{z) is 
related to the bulk pressure and the surface tension. We 
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FIG. 2: Normalized electric potential <l>(z)/A<l> and compo- 
sition (j>{^) (upper panel), and normalized ion densities ci{z) 
and 02(2:) (lower panel), where x = 3, ^ = 4, (71 = —g2 — 10, 
and cia = cip — 2 x 10~*. For this pair of hydrophilic and 
hydrophobic ions, a microphase separation forming a large 
electric double layer is apparent. 



shall see that the discontinuity of ui across the interface 
— i^f} is nearly equal to the minus of the pressure dis- 
continuity pp — Pa in equilibrium. For a planar interface 
we have lUq, = uj^ = lOoo ■ 

Using Eq.(2.1) and eliminating /i.^ with the aid of 
Eq.(2.8), we obtain 



""^^^E'-p^ (2.17) 



2 ■ ■ »7r 

where n = X^j"-* i^ ^he total ion density. Since lo{z) 
Woo as z ^ ±00 for a planar interface, h is expressed as 

h - [/o(0a) - /o(0/3)]/A0 - rAn/A<^, (2.18) 

where An = - Ufj. With Eqs.(2.13) and (2.18) the 
bulk volume fractions and 0^ are determined for given 
bulk ion densities (see Appendix A) . The surface tension 
is expressed as 7 = / dz[uj{z) — tOoo] (see Appendix B). 
It consists of two contributions as 7 = 71 + 70 with 

71= [ dz[fo{4>) + ^<l>'^ 'h4,-Tn-Ca], (2.19) 



7e = y dzp<i> = y dzeicP)E^ 



(2.20) 





-0.01 
-0.02 

-0.03 
-0.04 
-0.05 
-0.06 
-0.07 
-0.08 



\ ' ' ' ' 


p 




^ 2 


- \ 


-ra 








9,= 10 




o=-10 






Ay a^/ 













0.0004 



0.0008 



0.0012 
''la 



0.0016 



0.002 



FIG. 3: a^A'y/T and a^A'yi/T as functions of cia for strongly 
hydrophilic and hydrophobic ion pairs, where A7 — A71 — |7cj 
and A71 — 71—70. The latter is very close to —a^F in 
agreement with Eq.(2.22). In this case |7e| is considerably 
larger than |A7i|. The parameters are the same as in Fig. 2. 




FIG. 4: a^Aj/T (left) and a^F (right) vs 32 for cic = 2x 10"" 
(a) and 10""^ (b), where x = 3, A = 4, and gi = 10. Anions 
are hydrophilic for g2 > and hydrophobic for t;2 < 0. The 
Gibbs relation A7 = — TF holds for hydrophilic pairs, but 
does not hold for —32 ^ 1 due to increasing |7e|. 



where (j)' — d(j)/dz and Ca ~ foi4'a) — h(j)a — Tua- From 
Eq.(2.18) the integrand of Eq.(2.19) vanishes in the bulk 
a and (3 regions. The lower and upper bounds of the in- 
tegrations in Eqs.(2T9) and (2.20) are pushed to infinity. 

Here it is convenient to introduce the excess adsorption 
F of the ions onto the interface by 



F = dz\n — Ua — 



A0^^ V-ayj 
Zint rL 

dz{n — na)+ / dz{n — nij). (2.21) 

In the first line the integrand vanishes far from the inter- 
face. In the second line we suppose a finite system in the 
region < z < L with L 3> ^ and determine the interface 
position by Zint — Jq dz[(j){z) — (j)p]/A(j). In Appendix A 
we shall see that 71 is related to F at low ion densities by 



A71 =71-70 



-TF, 



(2.22) 
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where 70 is the surface tension without ions. The so- 
called Gibbs adsorption equation^^ relates the surface 
tension and the excess adsorption by {d'y /d\nn)T — 
—TV, where n is the number density of the doped par- 
ticles (the ion density in our case). This equation yields 
Eq.(2.22) if r oc ri and if 7c is neglected. In our problem, 
however, Eq.(2.22) can be derived even if T is not a linear 
function of n (as in Fig. 3). 

The 7o in Eq.(2.20) is the minus of the z integra- 
tion of the electrostatic energy density /d — eE^/Sn. 
It arises from the last two terms in Eq.(2.17) with the 
aid of / dzp^ = J dze{(t))E'^ / An . In all the examples in 
our previous work^^^, it was at most a few percents of 
A7 = 7 — 70- However, in the low density limit, we shall 
see that 70 is proportional to the square root of the ion 
density. Moreover, its magnitude will turn out to be en- 
larged for strongly hydrophilic and hydrophobic ion pairs. 
This square-root dependence can be understood easily 
as follows. For simplicity, if \gi\ and [52 1 are of order 
unity, the screening lengths in the two phases, are of the 
same order and ^{z) changes from $q to <I>^ smoothly 
on the scale of the screening length ~ kJ^. Thus 
E{z) ^ A<i>KQ for \z\ < l/Ka and 7e ^ -(Afi))^eQKQ,/47r. 



III. TWO SPECIES OF IONS 

In the presence of two species of ions, where Zi > and 
Z2 < 0, we have the common density ratio riia/nip = 
n2a/n2i3- In terms of giA.(j), the ion density ratios and 
the Galvani potential difference are given by 



"1/3 

T 



exp 
9i 



\Z2\g1 + Zig2 

Zl + \Z2\ 
92 



Acj) 



Z, 



-A0. 



(3.1) 
(3.2) 



The Galvani potential difference is created by an electric 
double layer at the interface. There is no electric field for 
the symmetric case g\ — g2 in our theory. The potential 
$(z) changes on the scale of in the a region and on 
the scale of in the /3 region far from the interface at 

low ion densities. When 00 (as in the water-air 

case), ^(z) — becomes very small around the interface 
changing slowly in the /3 regioni^, as has been assumed 
in the literatureii^i^. 



Numerical results 



the solvation parameters gi and the ion densities in the 
following figures. 

In Fig. 2, we show a set of equilibrium profiles for g\ — 
—g2 = 10 in the monovalent case'^^. Here c\a and c\p 
coincide from Eq.(3.1) and is set equal to 2 x 10"'*—. The 
interface thickness ^ is of order 5a, while $(2:) changes on 
the scale of k'^ ~ k~^^ ~ 10a. Here A7 ^ -O.G4ira-2 
and r — 0.014a^^. We can see a marked growth of the 
electric double layer, which is enhanced with increasing 
g\ for the case g\ = —g2- The ion density ci -I- C2 has 
a deep minimum at the interface position, for which see 
Fig. 5 also. In our previous worki^, we obtained milder 
ion profiles for gi = —g2 = 4 and x = 2.3. 

In Fig. 3, we examine how the surface tension 7 is 
decreased with increasing ciq,(= ci^), where gi = —172 = 
10 as in Fig. 2. We notice the following, (i) The changes 
A7 = 7 — 70 and A71 =71—70 are both proportional 



1 /2 J. z 

to at small cia- Here |A7|/c^^^ is of order unity, so 
A7 is appreciable even for very small ciq. On the other 
hand, for hydrophilic ion pairsi^, A7 = A\CiaT/a^ > 
with the coefficient A\ being of order unity, which is 
consistent with the well-known surface tension increase 
for water-air interfaces with salt, (ii) Comparing A71 
and A7 — A71 -l- 7o, we recognize that 70 is a dominant 
negative contribution in this case, (iii) We confirm that 
Eq.(2.22) holds excellently (see Appendix A). 

In Fig. 4, we display A7 and F as functions of 172 in the 
range [—10, 10] at gi = 10. For 32 > 0, we have the usual 
behavior A7 > and F < 0. For 172 < 0, their signs are 
reversed and their magnitudes are increased dramatically. 
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B. Analysis using the Poisson-Boltzmann equation 
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We numerically seek equilibrium interface solutions of 
Eqs.(2.8) and (2.10) for C = x and ei/sc = 0.8 as in 
our previous work^^. Hereafter cia — C2a and cip — 
C2/3 in the monovalent case. We set x = 3 and A — A, 
except in Fig.l. We have 70 = 0.497T/a^ without salt 
and A(j) = 0.93 at very small ion densities (see Appendix 
A) at X = 3. The dielectric constant Sa in the a phase 
is twice larger than that sp in the P phase. We vary 



FIG. 5: Normalized electrostatic energy density 
a'^£{(f>)E{z)'^ /8tvT and ion density deviation a"^(n(z) — Ua) 
calculated numerically, which are compared with the Poisson- 
Boltzmann solution a^{npB{z) — ria) in Eq.(3.4). The z 
integration and division by —a^jT of the first quantity is 
equal to 7c in Eq.(2.20) and that of the third quantity 
is equal to 7^^ in Eq.(3.7). Here x = 3, A = 4, and 

Clc = Cl;3 = 2 X 10"". 
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FIG. 6: Numerical a^^^jT in Eq.(2.20) and approximate 
a^'y^^ /T in Eq.(3.7) as functions of cia, which agree at small 
cia(< 10-*). 

Since 7e is dominant in Fig. 3, we seek its approx- 
imate expression to examine its overall behavior. To 
this end, we consider situations in which the image in- 
teraction is not crucial, the ion densities are very low, 
and the ions are monovalent. Then, far from the inter- 
face position \z — Zjntl > the ion distributions obey 
the nonlinear Poisson-Boltzmann equation. That is, for 
the normalized potential U{z) = e$(z)/r, we assume 
dU{z)'^ /dz'^ = K|sinh(t/(z) — Up) in the region z > Zi^t 
and dU{z)^/dz^ — n\ sinh([/(z) — Ua) in the region z < 
Zint, where Uk = Q^k/T and kk = {^■nriKe^ / ekTY^"^ 
[— Debye wave number) in the two phases, K — a and 
p. Here we are supposing the thin interface limit ^ — > 0. 

The resultant electric potential $pb is then given by 



2T 

^pb{z) = $k H In 

e 



1 - d^e"'^'^l^"^i"tl 



(3.3) 



where if = a in the region z < Zjnt and K — P in the re- 
gion z > Zint. The coefhcients da and djs are determined 
from the continuity of <i>pB and ed^pB/dz at z = Zjnt. 
In this approximation, the electrostatic energy density 
divided by T and the ion number density coincide as 



SttT 



= "-pb(z) -riK = 



[1 - d2^e-2KKl^-^i„tl]2 



(3.4) 



where E'pb — —d^pB/dz is the electric field and 
npBiz) — UK cosh[C/ (z) — Uk] is the number density with 
riK = niK + n2K being the bulk densities. The potential 
value at the interface position is given by 



T 



[4>a - $(zint)] = In 



1 + 5e^^*/2^ 



1 + 6e-«^*/2T 



(3.5) 



where b = ef^Kp/eaHa = npKa/riaKp with Eq, and 
being the dielectric constants in the bulk phases. In 
our scheme, Eqs.(3.1) and (3.2) give eA$/2T = (51 — 
g2)Acj)/4 and 



(3.6) 



Note that A$ and b are independent of the ion density. 

In Fig. 5, we compare numerically calculated n — 
Ua and eE^/Stt with the Poisson-Boltzmann solution 
eE^-g/SirT = npB - ria in Eq.(3.4). Here, ^ = 3, A = 4, 
and cia — cijs — 2 x 10^^, for which the image inter- 
action is not severe as shown in Fig. 1. Far from the 
interface the Poisson-Boltzmann solution is a good ap- 
proximation for these two quantities, while at z ~ Zmt it 
neglects the double-layer structure and cannot describe 
the deep minimum of n — n^. We notice that the z in- 
tegration of eEp^/SnT = npB — Ua is close to that of 
eE^/Sn but is a few times larger than that of F in this 
case. Hence we approximate 7e in Eq.(2.20) as the in- 
tegral of —eEp^/8nT = — rtpB + denoted by 7^^. 
Some calculations yield 



7, 



PB 



T 



2na 



1 + 6-^1 + 52 + 2& cosh 

As{naltBo.f'\ 



2T J 



(3.7) 



where b is defined by Eq.(3.6) and isa = e^/SaT is the 
Bjerrum length in the a phase. In the second line, the di- 
mensionless coefficient Ag is determined by b and eA$/r 

and is independent of the ion density ric, so 7^^ oc —n]/'^. 
Thus we obtain the result A7 cx — na^/^ at low ion densi- 
ties for asymmetric ion pairs gi ^92- For small eA$/T, 
we have 7c — (A$)2ec<KQ/47r in accord with the argu- 
ment at the end of Section 2. Also the excess adsorption 
F should exhibit this ion-density dependence, since the 
long-range tail of n — and that of npB — Ua should 
coincide. 

In the case gi > —g2 ^ 1, we have eA$/T = {gi + 
|52|)A0/2 and b < e'^^*/^^ in Eq.(3.7). Then, 



A.. 



r-1/2 



(e/3/£c 



\l/4g|g2|A0/4 



(3.8) 



which grows with increasing gi and \g2\- In Fig. 3, we 
notice that the coefficient in front of c\(^ in Aja^/T can 
be of order unity, which is equal to ^^(77/4^)^/^ from 
Eq.(3.7). In Fig. 6, we compare this approximate A7^^ 
and the numerical A70 for gi = —g2 = 10 as functions of 
cia- Agreement is excellent for cia ^ 10"'', while |A7^^| 
is larger than |A7e| for larger cia- 



C. Hydrophilic ion pairs 

On the basis of the Poisson-Boltzmann theory, we may 
also examine the usual case of hydrophilic ion pairs, 
where gi and g2 are both considerably larger than unity 
with gi > g2- (Here A7^^ = for 51 = (72.) Then we 
have Ha ^ np and 6 <C 1 not close to the critical point, 
leading to 



A7^^" ^ -2r[cosh(eA$/2T) - l]6n„/Kc 



(3.9) 



where bna/na = np/np oc np^^'^. Here the right hand 
side is the integration result in the /3 region, since that 
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in the a region is smaller by h. From Eq.(3.5) <I>q — 
$(z) c>c 6 in the a region. If (gi — 172)^^ > 1, we find 
|A7P|aVr - e-9^^'^/2(ci„/A)i/2, where A is defined 
by Eq.(2.6). Thus, together with positive A71, we obtain 
the following expression, 

A7 = -A.TK/^Ba)'/' + KTn^- (3.10) 

The coefficient A, is smah here (~ {epj ea)^l'^e-^'''^'^l'^). 
The second term is of the well-known form accounting for 
the ion depletion near the interface. With this expression, 
A7 should exhibit a small minimum given by 

(A7)„,in = -T\,n^ (3.11) 

at n = Urn = {AsjlXsY I l^a- As an example, let the 
ion concentration giving this minimum be 1 mM in the 
water-rich phase^i. Then we obtain yl^ = 1.2 x lO^'^ by 
setting As = SA and £bq = 7A. 

For water-air interfaces, Jones and Ray^ found a neg- 
ative minimum in A7 of order — 10~*7o. Their data can 
well be fitted to Eq.(3.10) with A^ ^ 10"^. However, 
we have assumed appreciable ion densities even in the 
less polar fi region. That is, our one-dimensional calcu- 
lations are justified only when the screening length 
in the /? region is much shorter than any characteristic 
lengths in experiments, which are the inverse curvature 
of the meniscus or the wavelength of capillary waves, for 
example. In the literatureii^ii^iSiiS ions are treated to be 
nonexistent in the air region, so in our scheme we do not 
still understand the Jones- Ray effect. 

Here it is worth noting that Nichols and Pratt^ the- 
oretically derived the square root dependence of A7 as 
the low-density asymptotic law when ions are appreciably 
soluble both in the two phases. They calculated the ion 
density deviations decaying on the scale of the screening 
length far from the interface, which vanish as — > as 
in our theory. 

IV. THREE SPECIES OF IONS 

Next we consider the ion distributions in the case 
of three ion species for water-oil interfaces. As in the 
experiment^, we assume monovalent ions with charges 
Q\ = e,Q2 = — e, and Q3 = for the three species. 
Namely, the first species consists of monovalent cations, 
while the second and third species are monovalent anions. 
The ion distributions are very complicated and numerical 
calculations are needed for their analysis. 

Before presenting numerical results, we first examine 
the potential difference A$, which is determined from 
Eqs.(2.12) and (2.13)^^'^*. For example, in terms of the 
ion densities in the a phase, it follows the equation for 
A$ in the form, 

ci^e^-^^*/^ = C2ae(9i-''^)'^'^ + C3„e(''i-''^)^^ (4.1) 

where ciq = C2a + c^a- This yields Eq.(3.2) for c^a — 
in the monovalent case. 
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FIG. 7: Normalized ion densities 01(2), C2(z), and 03(2) of 
three ion species in the monovalent case for C2a = 2 x 10~^ 
and C3/3 = 10~^. The species 1 and 2 are hydrophilic as 171 = 
g2 ~ 10, while the third species is hydrophobic as — —13. 
Then C2/3 = and cs^ ^ 1. 




z/a 



FIG. 8: Normalized ion densities ci{z), £2(2), and £3(2) in 
the presence of three ion species with gi = —10, 32 ~ 12, and 
33 = —13 on a semi- logarithmic scale, resembling to those in 
the experimenlryi'. The third species does not penetrate into 
the a region. 

However, in the experimen1>i^, the third species 
(TPB~) was strongly hydrophobic such that c^a 03^ 
was realized, while the second species (Br~) was hy- 
drophilic. Supposing such cases, let us assume (72 > 
and gs < and choose C2a and C3/3 as control parame- 
ters. If we set 

X = exp[eA$/r - (g, - 32) A0/2] , (4.2) 
Eq.(4.1) becomes a cubic equation, 

X - X-^ = 2R[l - x^e^3'~g2)/^<t>Y (4 3) 
where i?(oc C3/3) is defined by 

R = e(^'i+9^)^'^/2^3/3/2c2a. (4.4) 
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The right hand side of Eq.(4.3) arises in the presence of 
the third species. We may well assume xVs3^S2)A0 ^ i 
for large |(;3| » 1 to obtain X ^ R + + or 

|;A$ = i(.9i - 52) + ln(i? + (4.5) 

If (51 + 32)^0 3> 1, we readily reach the new regime 
R ^ I even for small 03,3, where X = 2R holds and 

That is, we find 

|;A$ ^ 31 A(/. + \n{c3p/c2a) (4.6) 

for e-(si+92)'^'^/2 < C3(i/c2c < e-(fi+93)A0/2, In tjiis 
case, Cia = C2q, Ci/3 = 03^3, C2/3 ^ e"(9i+f2)'^"^c^„/c3/3 < 

C3/3, and C3a = e(9i+S3)'^"*c2^/c2a < C2a. 

In Fig. 7, we display the ion distributions in the pres- 
ence of three ion species in the monovalent case with 
Qi = e, Q2 = — e, and Q3 = — e. Since the absolute 
values of gi are taken to be large, we can see steep and 
complex variations of the ion distributions around the in- 
terface. The first and second species are both hydrophilic 
but the third one is hydrophobic as gi = g2 = 10 and 
53 = -13. Here x = 3, eA$/T = 7.92, and 7 = 
0.446T/a2, while A$ = for 03^3 = since gi = 52- The 
peaks of ci and C3 are conspicuous. This is the case dis- 
cussed around Eq.(4.6), since X^e'^S3-g2)A^ ^ ^ 5 ^ 
and i? = 2.7 X 10^. 

In Fig. 8, the first and third species are hydropho- 
bic but the second species is hydrophilic as gi — —10, 
52 = 12, and 33 = -15, where x = 3.2, eA$/T ^ 7.01, 
and 7 = O.eOOT/a^ (with 70 = 0.620r/a2 at x = 3.2). 
Here ci and C2 exhibit sharp peaks. The ion distributions 
in this case can be compared with those in the experiment 
by Luo et aZ.~, so the curves are written on a semiloga- 
rithmic scale as in their paper. The adopted parameter 
values are inferred from their experimental data. 

V. SUMMARY 

Effects of ions in polar fluid mixtures are very complex 
because of the presence of the electrostatic, solvation, and 
image interactions. Our continuum theory takes account 
of these interactions, though it should be inaccurate on 
the angstrom scale. Our main results are as follows. In 
Fig. 2, we have illustrated the singular ion distributions 
around an interface when hydrophilic and hydrophobic 
ions coexist. We have given the general expression for the 
surface tension in electrolytes in Eqs.(2.19) and (2.20). 
In Fig. 3, the resultant surface tension change A7 is 
proportional to the square root of the ion density in the 
dilute limit. This dependence arises from the fact that 
the electrostatic field changes over the distance of the 
screening length far from the interface. When the image 
interaction is not severe, the electrostatic contribution 7e 
to the surface tension can be estimated by the Poisson- 
Boltzmann result in Eq.(3.7), as demonstrated in Fig. 



6. For hydrophilic ion pairs, we propose the expression 
(3.10) consisting of the negative square-root and positive 
linear terms. In the presence of three ion species, the ion 
distributions are very complex as in Figs. 7 and 8, where 
Fig. 8 corresponds to the experimental result^^. 

Using salts composed of hydrophilic and hydrophobic 
ions, a large decrease of 7 far from the critical poin t^^'^^ , 
ion distributions near an interface^^, and a mesophase 
near the critical poin1>i^ have already been reported. We 
propose more systematic surface-tension measurements 
of water-oil interfaces with dissimilar cations and anions. 

In this paper, we have neglected the clustering of ions, 
(i) For hydrophilic ions, this effect becomes conspicuous 
for densities larger than IM in aqueous solution o^°i^^'^^ . 
It is not clear how hydrophilic and hydrophobic ions can 
form clusters with increasing their densities, (ii) It is 
worth noting that ion aggregation was predicted to oc- 
cur near the critical point (even among ions of the same 
species) because of a long-range interaction mediated by 
the critical fluctuations (see Appendix A of our previ- 
ous paper—), (iii) For ionic surfactants, the amphiphilic 
interaction with water and oil needs to be included. Re- 
cently we have found that the adsorption of ionic sur- 
factant and counterions onto a water-oil interface is dra- 
matically enhanced when these two species constitute a 
hydrophilic and hydrophobic pair producing a large elec- 
tric double layer^**. 

In future, we should study the surface tension near 
the critical point in the presence of hydrophilic and hy- 
drophobic ions with \gi\ 3> 1. The mesophase formation 
realized for such ion pairs near the critical point should 
also be investigatedi^iJJiiii. Dynamics of phase ordering 
in ionic fluids with the solvation interaction has not yet 
been explored. Transient relaxations under applied elec- 
tric field could also be studied. 
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Appendix A: Derivation of the Gibbs relation 
for dilute ion densities 

We here derive Eq.(2.22) in the Ginzburg-Landau the- 
ory for small n. Let 00 (z) and ho be the composition 
profile and the chemical potential difference without ions. 
From Eq.(2.10) the deviation S(j) — <j) — (po obeys 

i 

to linear order, where /q ~ d'^fo/dcj)^. For the free energy 
density in Eq. (2. 2) we have = [0-i(l-(/))-^-2x]r/a3. 
Obviously, the left hand side of Eq.(Al) vanishes to linear 
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order if we multiply it by (f)' = d(f)/ dz and integrate over 
z. We can show that the right hand side also vanishes in 
the same procedure from / dz{eiE'^ (/)' /8tt + ^' p) = and 
/ dz{hi^(j>' + J2i f^imZfn'^) = where the latter follows 
from Eq.(4.24) of our previous paper^-'^ and the primes 
denote operating d/dz. This orthogonality to (f>' is the 
solubihty condition of Eq.(Al)i^. From Eqs.(2.13) and 
(2.18) the deviation Sh = h — ho is calculated as 



5h 



-TAn/A^- 



The bulk volume fractions are expanded as 

A7l\ 



A(j) J 



//o(0 



(A2) 



(A3) 



where K — a and (3. The cjfKa £^re the bulk values without 
ions. From Eq.(2.19) we may expand 71 with respect to 
64> — 4> — (j)Q up to second order as 



71 



70 



rr 



dz 



A'^2 



C(2)-a-^(</>(^)-0a) 



(A4) 



where C = /o'(0)(^0)V2 and AC = Ca - C/3- The 
terms linear in 54) vanish from the interface equation 
ho — /o('*/'o) ~ C(j)o without ions. The third and fourth 
terms are of second order in 64> and, if they are neglected, 
we obtain Eq.(2.22). However, the corrections grow near 
the critical point. 

Appendix B: Grand potential and surface 
tension in incompressible fluid mixtures 

In the text, we have neglected the deviation of the 
number density nt = ua + ub of the solvent from a ref- 
erence density assuming the common molecular size 
{a A — OLE = a). More generally, the system is character- 
ized by the solvent densities, 



UA = nt4>, hb = nt{l - 4>), 



(Bl) 



in addition to the ion densities ni and 712- We assume 
(j) = nA/{nA + ns), so (/> is also the molar fraction. The 
density deviation Snt = nt — is small. Hereafter we 
suppress the dependence on T. 

For nearly incompressible fluid mixtures, the depen- 
dence on Snt may be accounted for if we replace fo{4') 
Eq.(2.1) by the Helmholtz free energy density. 



/o(0, nt) = fo{(f)) +ao+ ai(0)(5nf 



02(0) 



iSntf, (B2) 



where oq is independent of the densities (being dependent 
only on T), but ai and 02 depend on cj) (and T). Far from 
the interface, where the image interaction vanishes and 
the fluid is homogeneous, the pressure p is given by 



P 



-fo - flo + airit + a2Snt +T^^'i 



(B3) 



up to first order in 6nt. The last term is the contribution 
from ions. The coefficient 02 is inversely proportional 



to the compressibility {dnt/dp)T4>/'nt and is large such 
that a2(5nt is appreciable, while a2{5ntY /2 is small in 
Eq.(B2). The chemical potentials of the two components 
HA = SF/SnA and — SF/SnB are expressed as 

fiA = {I - <j)){h + Ah)/nt + ai + a26nt, 

Hb = —4'{h + /S.h)/nt+ai+a25nt, (B4) 

where h is given in Eq.(2.10) and 



(B5) 



We find h -f A/i = nt{iJ,A — IJ-b)- The expressions for 
the ion chemical potentials /ii and fj,2 in Eq.(2.8) are un- 
changed. For nonvanishing dnt we define the generalized 
grand potential density by 

W = /o - UAPA - nB^J.B - «iMi - '^2^2 

^ uj - ntPB - <P^h+ {fo- fo), (B6) 

where uj is the grand potential density in Eq.(2.17) in the 
incompressible case and use has been made of nA^^A + 
nB^J■B = ntfiB + (f){h + Ah). In the bulk regions uj should 
be equal to the minus of the pressure p. We assume that 
the electric field vanishes far from the interface and there 
is no contribution of the Maxwell stress tensor. 

In particular, when the fluid is separated by a planar 
interface, the four chemical potentials pA, ^J■B, ^J■l, and /i2 
become homogeneous. Here iv tends to a constant —poo 
as z —^ ±00. In the right hand side of Eq.(B6), we may 
replace nt by and fo — fo by oq for small Snt . Then, as 



assumed in the text, uj tends to Uoo = —Poo + n^pB ~ oo 
as z — > ±00 and the surface tension is expressed as 



7 = y dz[uj{z) +pa, 

= j dz[uo{z) — ujo 



(B7) 



where the first line is general and the second line is its 
incompressible limit Snt — 0. 

We also consider an equilibrium spherical droplet with 
radius R in one of the two phases, say, the a phase. In 
this case, the total Helmholtz free energy F is minimized 
for given system volume V and the total particle num- 
bers, Na, Nb, and Ni = N2. The densities inside and 
outside the droplet are then determined such that the 
four chemical potentials become homogeneous (see Chap- 
ter 9.1 of the book by the present author'^'^). We further- 
more minimize F with respect to R in the presence of 
the surface free energy Att^R'^. This yields the famous 
pressure difference Ap = p^ — Pf3 = 27/i? (the Laplace 
law). Thus there arises a difference in the bulk values of 
Lu given by 



a+ujf3 = Ap = 27/i?. 



(B8) 



which holds in the limit Snt — s- 0. As i? ^ 00, we have 

'-^a = = Woo- 
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